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I. INTRODUCTION 

Since the formulation of the first singularity theorems it is generally conceded that singular spacetimes are of 
fundamental importance in general relativity. Geometrically, a singularity is defined via the notion of (geodesic) 
incompleteness, a viewpoint which fits in the singularity theorems of Hawking and Penrose (see, e.g., ||], Chap. 8), 
forcing us to regard a singularity as some kind of singular boundary point of spacetime. Recently, as an alternative, it 
has been suggested to describe (mild) singularities as internal points, where the field equations are satisfied in a weak 
(probably distributional) sense (cf. Q). General relativity as a physical theory is governed by particular physical 
equations; the focus of interest is the breakdown of physics which need not coincide with the breakdown of geometry. 

In the context of conical spacetimes algebras of generalized functions Q have been used to overcome the problem 
of simultaneously dealing with singular (i.e., distributional) metrics and the nonlinearities of general relativity ilJ- 
These techniques allow to assign to the cone metric a singular energy momentum tensor supported on a submanifold of 
codimension two, which, by a result of Geroch and Traschen Q, is not possible within classical (i.e., linear) distribution 
theory. 

The main focus of this work is a (nonlinear) distributional description of the Schwarzschild spacetime. Although the 
nature of the Schwarzschild singularity is much “worse” than the quasi-regular conical singularity, there are several 
distributional treatments in the literature (|, 0,|^,|g,|^), mainly motivated by the following considerations: the 
physical interpretation of the Schwarzschild metric is clear as long as we consider it merely as an exterior (vacuum) 
solution of an extended (sufficiently large) massive spherically symmetric body. Together with the interior solution it 
describes the entire spacetime. The concept of point particles—well understood in the context of linear field theories— 
suggests a mathematical idealization of the underlying physics: one would like to view the Schwarzschild solution as 
defined on the entire spacetime and regard it as generated by a point mass located at the origin and acting as the 
gravitational source. This of course amounts to the question of whether one can reasonably ascribe distributional 
curvature quantities to the Schwarzschild singularity at the origin. 

The emphasis of the present work lies on mathematical rigor. We derive the “physically expected” result for the 
distributional energy momentum tensor of the Schwarzschild geometry, i.e., Tq = {x), in a conceptually 
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satisfactory way. Additionally, we set up a unified language to comment on the respective merits of some of the 
approaches taken so far. In particular, we discuss questions of differentiable structure as well as smoothness and 
degeneracy problems of the regularized metrics, and present possible refinements and workarounds. These aims are 
accomplished using the framework of nonlinear generalized functions (Colombeau algebras) |^, ^ and, in particular, 
the geometric approach taken in ||^, . 


The paper is organized in the following way: in section II we discuss the conceptual as well as the mathematical 
prerequisites. In particular we comment on geometrical matters (differentiable structure, coordinate invariance) 
and recall the basic facts of nonlinear distributional geometry in the context of algebras of generalized functions. 
Moreover, we derive sensible regularizations of the singular functions to be used throughout the paper. Section HI 
is devoted to a first approach to the problem; a detailed discussion follows in section IV: we comment on problems 
and obstacles associated with the direct approach. Finally, in section ^ we present a new conceptually satisfactory 
method to derive the main result. Overly technical calculations are shifted to various appendices. In the final section 
VI we investigate the horizon and describe its distributional curvature. Using nonlinear distributional geometry and 
generalized functions it seems possible to show that the horizon singularity is only a coordinate singularity without 
leaving Schwarzschild coordinates. 


II. PREREQUISITES 

To begin with, let us have a look at the conceptually much simpler problem of point charges in Maxwell’s theory and 
consider the Coulomb solution i of an extended spherically symmetric body. In an idealized picture the charged body 
is reduced to a point charge, and this way of looking at the problem has proven to be very fruitful, mainly due to the 
following two reasons: first, the function ^ G C^(M^\{0}), since also in naturally gives rise to a distribution 

on Reinserting this distributional potential into the field equation we obtain Ai = — 47r^, which has the clear 
physical interpretation as the charge density of a point charge. Second, also in accordance with physical intuition, 
the situation may be interpreted in terms of the following sensible regularization scenario: consider a regularization 
of the “singular” potential by any sequence of (say smooth) functions converging weakly to Then, by virtue of 
linearity of the field equation, distribution theory guarantees that the corresponding sequence of charge densities will 
converge weakly to —AttS, i.e., the density of the point charge. 

The general relativistic case is much more involved. Consider the Schwarzschild metric inside the horizon: extending 
the spacetime to r = 0 we are confronted with several distinct problems. First—according to the standard picture of 
general relativity—no manifold structure is given at the singularity r = 0, since the field equations are meaningless 
there within the smooth category. In addition, the differentiable structure of the extended manifold cannot be 
uniquely determined from the differentiable structure of the original spacetime. This problem is dealt with in the 
relevant literature by fixing some differentiable structure by hand, most often the one induced by Cartesians associated 
to Schwarzschild coordinates. 

In analogy to the Maxwell case, we want to regard the metric as a distribution on the whole extended spacetime. 
Now, the second conceptual problem is due to the inherently nonlinear nature of general relativity: no distributional 
meaning can be given to the field equations, since it is not possible to calculate the curvature from a distributional 
metric. In the literature, this obstacle is circumvented by using various—more or less—ad-hoc regularization ap¬ 
proaches in order to calculate a regularized Ricci tensor within the smooth category. Eventually, its distributional 
limit is computed and—via the field equations—a distributional energy momentum tensor is obtained. This tensor 
may then be interpreted as distributional source of the Schwarzschild geometry |, 0 0 HI, HI. However, using 
ad-hoc regularizations we are confronted with the problem of regularization independence of the results which may 
not be suitably addressed within this setting. 

In this work, while arguing form a related point of view, we are going to use a different apparatus to deal with the 
nonlinearities of GR: the theory of algebras of generalized functions gives us the additional flexibility and power of a 
rigorous mathematical framework in which distributions may undergo nonlinear operations. In particular, following 
the procedure of Q, we will first model the distributionally extended Schwarzschild metric by a generalized metric 
obtained by a suitable (and general) regularization procedure. Then, after entering the generalized framework (cf. 
i HI ) we may calculate all the relevant curvature quantities from the generalized metric and subject it to the field 
equations. Note that within the generalized setting the field equations possess a well defined meaning. Finally, we 
may descend to the distributional level for the purpose of interpretation using the concept of association (see below). 

Note that, in general, a regularization procedure depends on the coordinate system it is performed in (for a 
diffeomorphism invariant notion of regularization using paths of mollifiers see ©HI)- However, since we had to 
fix the differentiable structure beforehand this is no further restriction. Actually, geometric considerations play an 
important role: as shown below, a brute force regularization attempt does not lead to a sensible description of the 
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problem at hand. Indeed, we shall see that a satisfactory treatment of the distributional Schwarzschild spacetime 
has to use the Kerr-Schild form of the Schwarzschild metric (which fixes both the differentiable structure and the 
coordinates); moreover, it must be retained during the whole regularization process. Note that this is in accordance 
with physical intuition since in the Kerr-Schild form the radial coordinate retains its spacial character near the 
singularity which of course is not the case in Schwarzschild coordinates. 

In the remainder of this section we are going to introduce some mathematical prerequisites. First, we are going 
to shortly recall generalized tensor analysis and generalized curvature (in Colombeau’s so-called special setting). For 
all further details we refer the reader to [Q Second, we explicitly calculate the regularization of the relevant 
components of the metric tensor to be used throughout the paper. 


Nonlinear distribntional geometry 

The basic idea of Colombeau’s theory of generalized functions |^, Q is regularization by sequences (nets) of smooth 
functions and the use of asymptotic estimates in terms of a regularization parameter s. Let (Me)££(oq] with € 
C°°(M) for all s (M a separable, smooth orientable Hausdorff manifold of dimension n). The algebra of generalized 
functions on M is defined as the quotient Q{M) := £m{M)/N{M) of the space £m{M) of sequences of moderate 
growth modulo the space N{M) of negligible sequences. More precisely the notions of moderateness resp. negligibility 
are defined by the following asymptotic estimates (X(M) denoting the space of smooth vector fields on M) 

£m{M) := {(it,),, : VK CC M, Vfc G Nq G N 

V^i,... G X{M) : sup 1% .. .Lj, Ueip)\ = 0{e~^)} 

pGK 

JV{M) := {(it,),: VK CC M, Vfc,9 G No 

V^i,... G X(M) : sup 1% .. .Lj, it,(p)| = 0(e'?))} 

p&K 

Elements of G{M) are denoted by it = cl[(it,),] = (ue)e +Af{M). With componentwise operations G{M) is a 
fine sheaf of differential algebras with respect to the Lie derivative defined by L^u := cl[(L^ite)£]. The spaces of 
moderate resp. negligible sequences and hence the algebra itself may be characterized locally, i.e., it G G{M) iff 
uo ipa & GilpaiVa)) for all charts {Va,tpa), where on the open set ipaiVa) C K" in the respective estimates Lie 
derivatives are replaced by partial derivatives. Smooth functions are embedded into G simply by the “constant” 
embedding a, i.e., a{f) := cl[(/),], hence C°°{M) is a faithful subalgebra of G{M). On open sets of R" compactly 
supported distributions are embedded into G via convolution with a mollifier p G iS(M") with unit integral satisfying 
/ p{x)x°‘dx = 0 for all |a| > 1; more precisely setting Pe{x) = {1 /p{x/s) we have l{w) := cl[(ic * Pe)e\- In case 
supp(ic) is not compact one uses a sheaf-theoretical construction. However, in the special case of the functions to be 
treated in the context of the Schwarzschild metric this will not be necessary (see below). From the explicit formula 
it is clear that (on open subsets of Euclidean space) embedding commutes with partial differentiation. On a general 
manifold, however, there is no canonical embedding of V available; a suitable replacement (cf. jlj]) is provided by 
physically motivated modeling and the use of the notion of association (see below). Inserting p G M into u G G{M) 
yields a well defined element of the ring of constants (also called generalized numbers) JC (corresponding to K = K 
resp. C), defined as the set of moderate nets of numbers ((r,), G with |r,| = 0{e~^) for some N) modulo 

negligible nets (|r,| = 0{£^) for each m). Einally, generalized functions on M are characterized by their generalized 
point values, i.e., by their values on points in M,, the space of equivalence classes of compactly supported nets 
(Pe), G with respect to the relation p, ^ p(. dhipe^p'e) = 0(e”’') for all m, where dh denotes the distance on 

M induced by any Riemannian metric. 

The tJ(M)-module of generalized sections in vector bundles—especially the space of generalized tensor fields 
Gg{M) —is defined along the same lines using analogous asymptotic estimates with respect to the norm induced 
by any Riemannian metric on the respective fibers. However, it is more convenient to use the following algebraic 
description of generalized tensor fields 


Gl[M) = G{M)®r:{M), ( 1 ) 

where T^{M) denotes the space of smooth tensor fields and the tensor product is taken over the module C°°{M). 
Hence generalized tensor fields are just given by classical ones with generalized coefficient functions. Many concepts 
of classical tensor analysis carry over to the generalized setting , in particular Lie derivatives with respect to both 
classical and generalized vector fields. Lie brackets, exterior algebra, etc. Moreover, generalized tensor fields may also 
be viewed as t/(M)-multilinear maps taking generalized vector and covector fields to generalized functions, i.e., as 
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0(M)-modules we have 


In particular a generalized metric is defined to be a symmetric, generalized (0, 2)-tensor field gab = cl[((/a{,e)e] (with 
its index independent of e and) whose determinant det{gab) is invertible in g{M). The latter condition is equivalent 
to the following notion called strictly nonzero on compact sets: for any representative det(5ab£)E of det(gah) we have 
Vif CC M 3m G N : infpgif |det(ga6£)| > e™ for all e small enough. This notion captures the intuitive idea of a 
generalized metric to be a sequence of classical metrics approaching a singular limit in the following sense: gab is a 
generalized metric iff (on every relatively compact open subset V of M) there exists a representative {gabs)s of gab such 
that for fixed e (small enough) gabs (resp. gabe\v) is a classical pseudo-Riemannian metric and det(5a6) is invertible in 
the algebra of generalized functions. A generalized metric induces a t/(M)-linear isomorphism from 5 q(M) to t/°(M) 
and the inverse metric := c\[{g~f^Ye] i® ^ well defined element of g^lM) (i.e., independent of the representative 
{gabe)e)- Also the generalized Levi-Civita connection as well as the generalized Riemann-, Ricci- and Einstein tensor 
of a generalized metric are defined simply by the usual coordinate formulae on the level of representatives. 

Finally, the setting introduced above displays maximal consistency (in the light of L. Schwartz impossibilty result 
[p^ ) with respect to smooth resp. distributional geometry most conveniently formalized in terms of the notion of 
association. A generalized function u G g{M) is called associated to zero, u « 0, if one (hence any) representative (us)e 
converges to zero weakly. (In a sloppy fashion we shall often write Ug ~ 0.) The equivalence relation u ~ v u—v « 0 
gives rise to a linear quotient of g that extends distributional equality. Moreover we call a distribution w G V'(M) 
the distributional shadow or macroscopic aspect of u and write u « ic if for all compactly supported n-forms iz and 
one (hence any) representative (us)s 


lim 


Ugly = {w,v), 


M 

where (, ) denotes the distributional action. By (|^) the concept of association extends to generalized tensor fields in 
a natural way. 


Regularizations of the singular functions occurring in the Schwarzschild metric 


The two most important singular functions we will work with throughout this paper (namely the singular compo¬ 
nents of the Schwarzschild metric) are ^ and {r = ||af|| ; c a positive constant). Since ^ G L;g^(M^), it gives rise 
to the regular distribution ^ G By convolution with a mollifier p (adapted to the symmetry of the spacetime, 

thus chosen radially symmetric) we embed it into the Colombeau algebra I/(K^) 


4) = (r*ft) =: (r)e. 


r r r 

Using radial symmetry of the convoluted function and inserting Ps{r) = ^p(|) we obtain 

rje 

A-tt r 

dt t p{t). 


( 2 ) 


(-)£ =— f dtY p{t) 
r r J 


47r 

£ 


( 3 ) 


r/e 


It is easy to confirm that (i)^ = a{^)s = ^ on R^\{0}, and at the origin we have (^)s = — /q°° dt tp{t) . 

r=0 

In contrast to the function is not in L)^^(R^). A canonical regularization (in the sense of Gelfand-Shilov 
[p^) is the principal value G I?'(]R^) which can be embedded into 5(R^). 


1 


I 


vp{ -) G !?'(»•') 


r — c r — c 

Making use of vp(-^) = -1; log |r - c| we obtain i(vp(-^))(a:) = 


i(vp(^—)) =: (vp(^—)), . 
r — c r — c 


( 4 ) 


d f 1 d f 1 

J ^^^ \x-y\ ~^\pAy )“ \x-y\ \^~y\ -^IpAv) 
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and finally for v > c 


I —C LXJ 

1 A-tt C A-tt C 

i.{vp{ - )){x) = - / dspe{s)s'^ -/ dsps{s)s'^ + 

r — c r — c J r J 

0 r — c 

0 

OO ^ 

/—0 

r-c 

For 0 < r < c the roles of r and c are interchanged and at the origin we obtain vp(^^) 


( 5 ) 


= -Ho(£). 


III. A FIRST APPROACH TO THE PROBLEM 


In this section we present a first approach to the “Schwarzschild point mass problem”, thereby essentially following 
earlier treatments in the literature (||, |^, |l^ |l^). However, we are going to use the language of nonlinear distri¬ 
butional geometry introduced above in order to obtain a unified view, which will enable us to carry out a detailed 
analysis of the previous approaches in the next section. 

In the usual Schwarzschild coordinates (t, r > 0, 0, (j)) the metric takes the form 


= h{r) dt‘^ — h{r) ^ dVL‘‘' with h{r) = —1 -|- 


2m 


( 6 ) 


Following the above discussion we consider the singular metric coefficient h{r) as an element of C I?'( 

and embed it into by convolution with a mollifier. Note that, accordingly, we have fixed the differentiable 

structure of the manifold: the Cartesian coordinates associated with the spherical Schwarzschild coordinates in (|^) 
are extended through the origin. We have 


h{r) = -1 H- 

r 


L{h{r)) = he{r) = -I -|- 2m(-)£ G 




( 7 ) 


where is given by i)- Inserting into §) we obtain a generalized object modeling the singular Schwarzschild 
metric, i.e., 

ds^ = he{r) df - h,{r)-^ dr^ + dn^ . ( 8 ) 

The generalized Ricci tensor may now be calculated componentwise using the classical formulae 


K). = ( 7 ?\)e = \(^K + ^K'^=lAh, 


{R\)e = {R\)e = ^ 


1 + hr 


Due to the linear structure of Rg it is evident that it is associated to a constant times the d-distribution, i.e., 

1 


(Rgje = - Ahr = m A{-)r —AirmS (e —> 0). 


( 9 ) 

( 10 ) 

( 11 ) 


Investigating the weak limit of the angular components of the Ricci tensor (using the abbreviation $(r) = 
f siiidddd<f> we get (cf. appendix]^ 

[{R\)r^ d^x = ({rh'r -I- 1 -I- /ie)$(r)dr = 


y Snm 


r/s 


tp{t)dt] ^{r)dr = Sirm / dx^{ex) / tp{t)dt 


327r^m$(0) J dx J tp{t)dt — ^ 87rm(d|$) (e ^ 0). 
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Hence, the Ricci tensor and the curvature scalar R are of i5-type, i.e., 


i?g = R\ « — 47rm5 i?2 = ^3 ~ SirmS R « mnS. (12) 

Equations ( [l^ are obviously given in spherical coordinates. Strictly speaking this is not sensible, because the basis 
fields are not globally defined. Representing distributions concentrated at the origin requires a basis 

regular at the origin. Transforming the results for (i?j)e (i-e., (^) and ®) into Cartesian coordinates associated 
with the spherical ones (i.e., {r, 0, (j)} {x*}) we obtain, e.g., for the Einstein tensor 


G) « -8TTm6SlSl (13) 

Note that the use of the particular regularization (|^) is not essential here. We could have replaced (J^ by any other 
smooth ad-hoc regularization of h{r), as has been done, e.g., in [^, by setting h^^r) = —1 -I- -^^== 5 . Indeed, we can 
show that the results ( [l^ ) hold for all regularizations, i.e., for all sequences of the form he{r) = —1 -|- 2mse{r) h{r) 
(i.e., V Se smooth, Sg in V'). Eor the (0,0)- and (1, l)-components of the Ricci tensor the result follows from the 
special form of ( 0 ). Eor the angular components (cf. (10)) we write 


00 00 

2m [ r^( '^'^(^) -I- ) l>(r) dr = —2m [ dr - l>'(r) ^ 87rm$(0), (14) 

J r Jr 

0 0 


where in the last step we used the fact that ^ 4)' G I?(IR). 


IV. COMMENTS AND PROBLEMS 

In order to be able to calculate the curvature from the metric we must keep the regularization /ie(r) smooth on the 
entire spacetime. This fact—although somewhat hidden because we worked with spherical coordinates—is essential 
from the conceptual point of view. In fact, choosing a regularization hg{r) which is smooth only on R^\{0} is not 
sufficient to derive the result as is explicitly shown by the following counterexample. Set = — I-|-2toSe and ^ 
(with Oe —> I weakly) consisting of regular distributions, so that Sg G L;^^^(R^) with ^ e Moreover, we 

may require Oe(r) to be smooth on R^\{0}. Summing over (^) and ( |l^ we get R^ = 2m(^ o" -I- ^ o(.). Choosing 
Oe(r) = (I -|- c [r® — 1]) we obtain for R^ different weak limits as the constant c varies, i.e., 

—> SirmcS. (15) 

For Oe = r~^ the situation is even worse. Although Oe G V e and Oe ^ 1 G P' as £ ^ 0, i?E does not converge 
weakly, so that we obtain no distributional result whatsoever. 

Nonetheless, similar non-smooth regularizations have been considered in the literature. In these cases the desired 
result (p^ can only be produced by means of a clever choice of explicit formulae; in particular, Oe = in |]^ and 
Oe{r) = 0(r — e) in [|2). The authors of ||] have shown that the result ( |l^ may be reproduced as long as Oe\r=o = 0. 
However, the conceptual problem remains untouched: Re[he] can only be derived for smooth regularizations 
distributions cannot be used as an input for nonlinear operations. 

Prior to a more detailed investigation of the choice of regularization, we briefly comment on two more attempts 
in the literature. In a regularization of the metric using thin shell solutions is investigated. The limit (e ^ 0) 
corresponds to a shrinking of the shell. However, the shells can only be placed outside the horizon (of a black hole 
with identical mass). This implies that a shrinking of the shell must be coupled to a decrease in mass: m converges 
to zero in the limiting process, so the obtained results should either be considered trivial {R ^ mS \m=o = 0) or be 
rejected completely. 

In fill the authors claim to have found different results for the curvature quantities by regularizing the Schwarzschild 
metric in a different coordinate system. They study the interrelations of regularizations and coordinate transformations 
for this particular problem. However, some details are not overly convincing. If we choose a new radial coordinate 
f such that r = Af -I- a with a = 2m, then f = 0 does not describe the Schwarzschild singularity. Instead, r = 0 
corresponds to the coordinate singularity at the horizon r = 2m, but shrunk to one point. Obviously, we should not 
compare the outcome of these considerations with our former results. 

We now begin with an in-depth analysis of certain aspects of the regularization procedure commencing with the issue 
of componentwise regularization and invertibility of the regularized metric. According to (^ in section H, regularizing 
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a tensor such as the Schwarzschild metric comes up to regularizing each distribution-valued component separately. 
Following this rule we obtain a regularized metric slightly different from (^, namely 

dSg = hs{r) dt^ — {h~^)e{r) dr"^ + r'^ dU^ . (16) 

Since cl[/ie]cl[(/i“^)e] 1 g the determinant of the regularized metric is no longer identically one. (This, 

in fact, does not come as a surprise; cf. Schwartz’ impossibility result [^.) However, the product is preserved in 
the sense of association, i.e., he{h~^)e ~ 1. Analogous issues arise from the inverse metric: embedding also g~^ 
componentwise into Q we obtain regularized objects, and (g~^)e, which are only inverse to each other in the sense 
of association. Taking a different viewpoint, however, it is comparatively easy to avoid these issues: on the classical 
level the Schwarzschild geometry is uniquely determined by the set of variables {gtt, grr, gee, gipip}, or, e.g., equivalently 
by { 9 tt, 9 ee, 9 (p(p,detg}. Embedding the second set of variables into G leads directly to the regularization (||) used 
above; no invertibility problems arise at all since det g^ is forced to equal one. 

Finally we return to discussing the problem of smoothness of the regularized metric from a different, more geomet¬ 
rical point of view. We regard this problem to be so essential that in the next section we propose an approach entirely 
different form the one taken so far. 

In fact, the regularizations used so far (as all the other regularizations in the relevant literature) do not provide 
a smooth regularized metric tensor. This fact is hidden again by the use of spherical coordinates. In Cartesian 
coordinates pertaining to (r, 0, </>)—which we used to fix the differentiable structure of the extended manifold at 
r = 0—however, it can be explicitly seen from the form of the metric 

ds^ = h{r) dt^ + ddP' — (1 + dx'^dx^ . (17) 

In order to obtain a smooth regularization it is not sufficient to merely regularize h{r). In fact, we must embed the 
singular coefficient (1 -I- /i(r)“^)^^ as a whole into G- Apart from technical difficulties this approach should provide 
a smooth regularized metric ds^. However, we have reached an impasse: the regularized metric will not be invertible 
at some distinct value rg of the radial coordinate, where rg ^ 0 (e ^ 0). This will be shown in the remainder of this 
section. 

As shown in appendix ^ the regularization of dH) takes the form 

dsl = hg{r) dt^ + (1 — eLs{r)) dr^ + (1 — bs(r)) r'^dVi^, (18) 

with a£(0) —> I (e ^ 0). This implies that the rr-component of the regularized metric (p^ ) is positive at r = 0 (at 
least for small e), because (grr)e(O) = (1 — a£(0)) ^ | (e ^ 0). On the other hand, (1 — ag) approximates —h~^, 

i.e., {grr ^c{x ^ 0) ^ ' 2 m _ r ^ b {^ ^ b). So w^e couclude that at some value cg of v the smooth function {grr ^€{^) 

must have a zero at least for small e. (Interestingly enough, this is not the case for negative masses since — 2 m-r- 
positive then). In other words, this means that the regularization of the metric © degenerates at some radius rg. 
Evidently, rg ^ 0 as e ^ 0. 

Note that the occurrence of this radius of degeneracy is neither due to the fact that we choose the particular 
regularization (p^, nor is it possible to avoid it by giving up spherical symmetry. To see this in some more detail 
consider the spatial part of (© (set h{r) := h ^(r)) and consider a certain class of regularizations 

dSg = dx"^ — (1 + h^{r)) dx''dx\ (19) 

where he{r) denotes an arbitrary regularization of h{r). However, for ds^ to become smooth, we must require that 
hg{r) be — l + 0{r'^) for (r —> 0). Now, an arbitrary regularization of not necessarily respecting spherical symmetry 
is obtained by adding zero-sequences to (p^. 

= ddP — (1 + he{r)) dx^dx^ + (ay)£(a:) dx^'dx^. (20) 

For special cases of {aij)e it is easy to show that ( po|) is degenerate. Choose, e.g., {aij)e such that only (ai 2 )£ =: he{x) is 
non-vanishing. The determinant of ( p0| ) at x = (0,0, r) is equal to —hs{r){\ — he{xY). As —^£(0) = 1 and —hs{r) < 0 
for small £ and finite r, there exists a radius (with r£ ^ 0 for e ^ 0), such that det(( 7 £) = 0. Again, we observe 
degeneracy. 







V. THE KERR-SCHILD APPROACH 


To begin with let us summarize what we have done so far: we considered regularizations of the Schwarzschild metric 
(using the language of algebras of generalized functions) to calculate the (distributional) curvature at the singularity. 
The regularizations used were essentially based on Cartesian coordinates associated with the spherical Schwarzschild 
coordinates. However, it turned out that all these regularizations were either non-smooth or not invertible. Smoothness 
and invertibility mutually exclude each other in this context. Hence, we are going to take another more geometrical 
view-point in this section. The main idea—following |l^ —is to use the Kerr-Schild form of the Schwarzschild 
metric. Retaining this preferred structure also during the whole regularization process will enable us to derive the 
physically desired result in a rigorous manner. 

A metric belongs to the so-called Kerr-Schild class if it can be written as 

gij = gii + fkikj with k^ki = 0. (21) 

Here, the null vector field k G X is normalized (fc*^ = 1) and / is a smooth function. Exploiting the Kerr-Schild form, 
some curvature quantities take a particularly simple form, e.g., 

R = dMk‘^k^)- ( 22 ) 

The Schwarzschild metric is a member of the Kerr-Schild class. In fact, transformation to Eddington-Finkelstein 
coordinates {t = t + 2m log |2m — r|) yields 


= —dP + dr"^ + r'^d^P ‘^^{dt — dr^. 

r 


Evidently, ( p^ ) is of Kerr-Schild form, gij = rjij + fpkj , with 


/ = —. 
r 


(23) 


(24) 


Analogously to section HI, we regard / and k as distributions on R'^. By this we again implicitly fix the differentiable 
structure: the coordinates {t,x'') are extended through the origin. 

We now proceed by regularizing both f and k. Indeed, this is necessary due to the fact that not only the profile 
function / is singular, but also the null vector field k is non-smooth. Recall that, on account of the nonlinearities in 
R[gr], (pH) can only be derived for smooth functions; it is inaccessible for distributional input. (Note the analogy of this 
situation with the one encountered in (H)). Hence, f and k are chosen to be the fundamental va riab les characterizing 
the metric (compare with the remarks in Section^). Regularizing the function / as in Section HI gives 


, N 2m 

/(r) = — 


Pf{r)) = fe{r) = 2m(i) 


The regularization of ki is carried out in detail in appendix yielding 

h(r) = ^ ^ i(^) = h, = (^), = X. F,(r) (z = 1,..., 3), 


(25) 


(26) 


where is given by (C3). Note that, for the moment, fco = I is embedded trivially into Q. Collecting the results of 
(p5|) and d^) we get the regularized metric 


ds1 = (—I -I- fe)dP — 2rfgFg dtdr -b (1 -f fe)dr^ + r^d^^ 


(27) 


where, for simplicity, t has been replaced by t again, and abbreviates r'^f^F^. Unfortunately, (|2^) is no longer 
of Kerr-Schild form. This can be shown indirectly: assuming that (E^ is Kerr-Schild, i.e., assuming that ds^ = 


-dt'^ 

dt \2 


- dr^ 


' /g {dt' — drY can be achieved by a transformation t t'{t, r ), it follows that /' = 


. As a consequence, the denominator 1 — fe + fe must be a strictly positive function. In fact, in the 


( d ±,\2 _ 1 

\dt') - 1 -/,+/, ^ ^ 

sense of association, it is equal to one. However, at the origin r = 0, we obtain 1 — /g -|- /g = I — tp{t)dt , 

which is negative for small e as long as JY° tp{t)dt > 0, a contradiction. 

The fact that the embedding ( R^ is no longer of Kerr-Schild form bears a strong relation to the fact that smooth 
regularizations degenerate at a certain value of r (see Section IV): the determinant of ( p^ contains the factor 
1 — /g -b /e , which was shown to possess a zero. 
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Additionally, in analogy to the statements made at the end of Section IV, we note that the loss of the Kerr-Schild 
form does not stem from choosing to regularize the singular coefficient functions via convolution. On the contrary, it 
may be shown that any regularization of the metric displays this behavior as long as only the spatial components of 
k are taken into account. 

We will now take the announced geometrical view-point: we consider regularizations retaining the Kerr-Schild 
decomposition. This requires, in particular, that the regularized vector ke is still null. Thus, we consider the 
regularization 


kl=rF,k^ (t = 0... 3). (28) 

While the spatial components of coincide with (|2^), is only associated to 1 {k° = rF^ fv 1 ). As required, fcg 
satisfies the condition klkei = 0. Note that, in order to obtain the functions = and k ■ k are 

chosen as fundamental variables determining the geometric structure of the spacetime. 

Using ( p5| ) and (^) the regularized metric takes the form 

9ij 6 ~ Vij “f ~ Vij “f (T fs) kikj = TJj^j k^kj. (^1^) 

Obviously, is of Kerr-Schild form. Finally we have arrived at a regularization of the Schwarzschild metric which 
is both smooth and invertible (the inverse being — f^Fk^). This allows us to fully exploit the Kerr-Schild form, 
i.e., to use (p^), to obtain 

Re = dadb{fek‘^k’’) = ^ fe + fe + ^ fe ■ (30) 

To complete our program we calculate the weak limit of R^. The technically involved calculations are deferred to 
appendix p[ Finally we derive 

{Re)e « 8nmS. (31) 


The Ricci tensor can be treated in complete analogy to obtain the Einstein tensor 


(G“,Je « -87Tm6SS6°. 


(32) 


VI. THE HORIZON 


In this last section we leave the neighborhood of the singularity at the origin and turn to the singularity at 
the horizon. The question we are aiming at is the following: using distributional geometry (thus without leaving 
Schwarzschild coordinates), is it possible to show that the horizon singularity of the Schwarzschild metric is merely a 
coordinate singularity? In order to investigate this issue we calculate the distributional curvature at the horizon (in 
Schwarzschild coordinates). 

Examining the Schwarzschild metric in a neighborhood of the horizon, we see that, whereas h(r) is smooth, 
h~^{r) is not even (note that the origin is now always excluded from our considerations; the space we are working 
on is M^\{0}). Thus, regularizing the Schwarzschild metric amounts to embedding h~^ into Q (as done in (H)). 

dsl = h{r) dt^ - {h-^)e{r) dr"^ F (33) 


with h{r) = — 1 -I- 


2 m 


and {h ^)e(r) = —1 — 2 m[vp( 


1 


r — 2m 


)]. 


(34) 


Obviously, ( p^ is degenerate at r = 2m, because h{r) is zero at the horizon. However, this does not come as a 
surprise. Both h{r) and h~^(r) are positive outside of the black hole and negative in the interior. As a consequence 
any (smooth) regularization of /i(r) (or h~^) must pass through zero somewhere and, additionally, this zero must 
converge to r = 2m as the regularization parameter goes to zero (note the analogy to the situation in section IV). 

Due to the degeneracy of (^^, the Levi-Civita connection is not available. Consider, therefore, the following 
connection r(,^- G Q: 


^kj ~ [''i9)mk,j + l'{g)mj,k k{g)kj,m] 


Clearly, F coincides with the Levi-Civita connection on ]R^\{r = 2m, r = 0}, as i{g) = g and i{g ^) = g ^ there. 


( 35 ) 
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Unfortunately, F does not respect the regularized metric b{g) (^, i.e., b{g)ij-,k 0, e.g., t(g)oo;i = (1 — h{h~^)g)h'. 
However, compatibility with the metric L{g) is a priori ruled out by the following statement: there exists no connection 
whatsoever under which t{g) would be a parallel tensor. To show this, just look at (T denoting a not necessarily 
torsion-free connection) i(g)oo;i = '•( 5 ) 00,1 ~ 2L%i(g)oo ■ At the horizon 1 ( 5)00 = 0, so that 1 ( 5 ) 00 ;! |r= 2 m = h'(2m) = 
^ sense of association, however, the connection ( ^^ is in fact metric compatible: b{g)ij;k ~ 0. 

We now investigate the curvature pertaining to the connection (^), picking out i?oo£ as a characteristic example. 
The result of the calculations displays the following structure 


9 ^ 4 4 9 

^ m m .TU\ ^ N 

Rooeir) = vp^(- - + - 4—) + vp^{2—- 4—) + { - 3 -- 2—) 


C50 00 ^ 

= ^p'e{r)^cy +vp,{r)^dix^ - — 


/=2 


1^1 


8 m^ 


CX) 


1 


(36) 

(37) 


Here, the abbreviations vp^ = [vp{:p^^)]s and x = are used; ci, di and e; are constants. Equation holds 

for |x| < 1; the infinite sums converge in this case. 

If the horizon is excluded, Rme = 0 (mod N), because ( p5| ) coincides with the Schwarzschild Levi-Civita connection 
there. In the neighborhood of r = 2m we aim at comparing Rooe{r) with a Colombeau object of the type 
(/ a Schwartz function). To this end we choose a fundamental sequence = 2m -I- e'^oo converging to r = 2m and 
examine Roosirg) (use (|^ together with (||)). 

• q> 1: RoQe{re) = const -|- o(e‘?“^). 

• q < 1: Using (^) we find that vp^{rg) = ^ and vp'^{rg) = — _ 2 rra)^ (^’^ sense of generalized numbers). 

Inserting these results into equation we obtain Rooe{re) = 0. 

• q = 1: RoQg{rg) = const-|-o(l). 

Thus, Roog{rg) has the same asymptotic behavior as a sequence of the type /( ''^~^"^ ) (as £ ^ 0). As a consequence, 
the weak limit of Rooe{r) can be calculated easily, simply by evaluating f dr r^$(r)/( ’'~^™ ). Evidently, this expression 
vanishes as e ^ 0. Since analogous results hold for the other components of the Ricci tensor, we are finally able to 
state 


Rrje-0. (38) 

In other words: viewed as a distribution, Rij = 0 on R^\{0}, i.e., including the horizon. If we were courageous 
enough we could take this as a proof that the metric singularity at the horizon is only a coordinate singularity. 

We conclude this section with a remark on the connection (|^). Due to the degeneracy of any regularization of 
the metric (e.g. (|^) no canonical (Levi-Civita) connection could be defined. The choice of connection (|^ ) bears 
a strong relation to the regularized metric; however, there seems no way of telling if this choice is canonical in 
some sense and thus preferred to other choices. Despite this open question, at least it is clear that the connection 
(^) is a regularization of the Schwarzschild connection. Indeed, we could change our viewpoint: we consider the 
Schwarzschild connection (forgetting where it came from, i.e., forgetting about the metric), regularize its distribution¬ 
valued components and calculate the distributional curvature from it. Proceeding in this manner, we obtain the 
result (|^, i.e., the spacetime is weakly Ricci-flat (the origin was excluded from our considerations). 


APPENDIX A: MOLLIFIER INTEGRALS 


Throughout this paper we work invariably with radially symmetric mollifiers p(r) (cf. section H). Most importantly, 
we have the properties 


00 00 

J dtt^p{t) = ^ J dtt^^p{t) =0 (fc > 1). 


0 0 

We investigate multiple integrals involving the mollifier p(r) and powers of r: 

j,n+l foo 


nOO nOO ^n+1 /*00 -1 /•OO 

/ dtt" / sp{s)ds = -- / tp{t)dt-\ -- / dtt^^'^pU) 

Jx Jt n + l n+l 


(Al) 


(A2) 
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( [A2| ) holds for (n, k ^ —1), it is proven by simply per forming integration by parts. 
One of the most interesting cases resulting from (^) is n = 0 and x ^ 0: 


POO pco ^ 

/ dt sp(s)ds = — 

Jo Jt 47r 


(A3) 


APPENDIX B: EMBEDDING OF THE CARTESIAN COMPONENTS 


Referring to section IV we investigate 
' 1 + h{r)~^ 


XiXj\dx'^dx^ = \2m / i{\\^\) ZiZj pe{\\z + x\\) d'^ z\ dx^dx^ 


where f(q) = 


1 1 


2m —q 


. In order to simplify (Bl) we show the following relation 


1 + h{r)~ 


XiXj ) = Xi Xj Ce{x) for i j (cg smooth). 


Proof: Since both /(||i]|) and pe are even functions in Zi, we observe that 

XiXj) =2m f f{\\^\)ziZjPe{...,Zi,...)d^z = 0. 
0Ci=O J 


(Bl) 


(B2) 


(B3) 


We can conclude that 


— XiXj) = Xi c(.(x) 


(c(. smooth) 


(B4) 


Also, — XiXj)\x^=o = 0 , from which follows that c'^{x) = XjC^{x), yielding (^). Note, however, that the 

smooth function Cg{x) in (^) is not equal to —). ■ 

In the case i = j, equation (^) is no longer valid. We are able to calculate the m— component in the limiting case 
e ^ 0, i.e.. 


1 + h{r)~ 


x—0 


= 2m / f{\\^\)zfpei\\z\\)d^z -> - (e^O). 


(B5) 


Proof: Clearly, 2m J /(||z||)zfpg(||i)|)d^z is independent of the choice of the index i, so that we may substitute it by 
^ //(||z)|)||z|pp£(||.?||)(i^z. Obviously, this converges to ^ as e goes to zero. ■ 

The regularized metric ( |B1| ) is radially symmetric, R*g£ = {R a rotation), as long as radially symmetric mollifiers 
are used. Thus, it must be of the form of a general radially symmetric metric ds^ = o(r) dr^ + r^b{r) dfi^, hence 


1 + h{r) ^ 


XiXj 1 dx^dx^ = (tte — hs){r) — —dx'dx^ + bJr) dx 

I J, J. 


(B6) 


For the general radially symmetric metric (|B6| ) to be smooth {a£ — b£)(r) = O(r^) is required. We observe consistency 
with (B2) and conclude 


ag(r) = 6g(r) + 


(B7) 


At the ori gin r = 0 only the second term 6(r)remains relevant, since (a—6)(r)^^|r=o = C£XiXjlr=o = 0. Turning 
equation ( |B5| ) to good account we obtain 

6g(0)^i (e^O), ag(0)->i (e ^ 0). (B8) 

Combining with (p^), we finally obtain 

ds1 = ^^(r) dt^ + (1 — ae{r))dr‘^ + (1 ~ b£{r))r^did^ . (B9) 
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APPENDIX C: EMBEDDING OF k' 

We explicitly embed the radially outward pointing unit vector field ki = ^ (* = 1 • • • 3) into the Colombeau algebra, 
i.e., 


6(^)= 

r J \\x-x'\\ J PI 


\x + z\\ 


(Cl) 


Equation is of an analogous form as ( pl| ) in appendix]^. We may conclude that t(^) is a radially symmetric 
vector field. Moreover, despite (.(^) ^ XiL {^), we must still have (repeating (B2)ff.) 


,{^) = x,F,{x) (i = 1...3). 


(C2) 


Here, F^{x) is a smooth function and moreover, because this function must be radially symmetric, Fg(x) = Fg{r). 
This fact makes it possible to calculate i( —) explicitly. Take x = (0,0, r) and investigate t{ — ) = x^Fg[r) = rF^{r): 


rFe{r) = 


ykPe{r')d^x' 


r — r' cos O' 


= 2tt J r'^dr' ps{r') J d{cos0') 

= J 0(k - r'\ + (r + r')) + ^^^(k “ - (r + r'f)^ 


v/j .2 _|_ J.' 2 _ 2rr' cos 6' 


47r 


7 CXJ 7 > 

J dss"^ Peis) + ^ Jdsspeis)-^^ J dsS^Peis) 


(C3) 


Clearly, Fe{r) = ^ on K^\{0} and, moreover, Fe(r) i on the whole space. We can write the latter also in the 
form rFe(r) « 1. 


APPENDIX D: WEAK LIMITS FOR THE KERR-SCHILD CASE 


We investigate the distributional limit of (^). Inserting for fe, (^) becomes 


Rs = 


1287r^m / , 4 2,,, 4 


.[l]44],^-^[4],[2]^ 


(i6[1]3 + -[2],[4]^ 

+ ^[2],[1]2 + - 3pe[2]^ - 4rp,[l],[2], + 

- ^,PS]1 + ‘2^Peme[i]e + ^PeilUMe ) , 


(Dl) 


where [l]e := / dsspeis), [2]e := f dss^pe(s), [4]e := / dss'^pe(s). 

r 0 0 _ 

In order to compute the weak limit of (|Dl|), expressions of the form (^) and ( |D3|) b elow have to be investigated. 
(Note that the negative powers of r are compensated by the integrals so that both (p^) and (D3) are well-defined as 
r —> 0). 

r-4+2^+3v,(r) [lYe[2]i[4\i 


-9+2j+3i 


[1]; [2]/ [4]e 


ii+j + k = 2) 
(* + j + A = 3) 


(D2) 

(D3) 


Terms of the forms (]D2D and (D3) possess related distributional limits: 




-9+2i+3i 


[l]eM2]i[4]^), 


-k 


2j -k 3i - 6 


(r-^+2^+3V,(r)[l]rM2]i[4]^),+ 


(r-4+2^+3V,(r)[l]p2]p4]pi)^ . 


2j -k 3z — 6 
k 

2j -k 3z — 6 
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In order to show this, consider the distributional action on a test function $(£), i.e., f dr $(r)r 7+2j+3z pji |■ 4 J^fe^ 

_ 0 

Here, we have introduced $(r) := f smOd0d(/><P(x). Integrating by parts and using $'(0) = 0, the claim is established. 
Taking (j = 2; fc = 1) as an example, we obtain the following weak limit: 

y dr$(r)^ [ 2 ]^ 2 [ 4 ]^ ^ ^(0) J dxp(x)[2][4] + ^^(0) J dxx^p(x)[2]^ 

® $(0)|dxp(x)[2][4] + i^|.(0). 

OO CC X 

Here, [1] := / dttp{t) [2] := J dtt^p{t) [4] := f dtt‘^p{t) . 

X 0 0 

Eventually, we obtain the result (|^) for the distributional limit of (ill)- 
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